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Abstract 

We discuss the formation of graded morphogen profiles in a ceh layer by nonlinear transport 
phenomena, important for patterning developing organisms. We focus on a process termed tran- 
scytosis, where morphogen transport results from binding of ligands to receptors on the cell surface, 
incorporation into the cell and subsequent externalization. Starting from a microscopic model, we 
derive effective transport equations. We show that, in contrast to morphogen transport by extra- 
cellular diffusion, transcytosis leads to robust ligand profiles which are insensitive to the rate of 
ligand production. 
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An essential feature during the development of an organism is the emergence of different 
cell types. This process, called cell differentiation, is intimately linked to the position of the 
cells in the organism. About fifty years ago, Turing suggested the existence of molecules 
that provide positional information for differentiating cells Turing proposed that these 
morphogens self-organize into patterns which structure the future organism. Pattern for- 
mation based on reaction diffusion systems has been studied extensively Today, 
a number of proteins have been identified to act as morphogens. In contrast to Turing's 
original suggestion, morphogens are secreted from a spatially localized source and form a 
stable, long-ranged concentration gradient in the adjacent tissue j^. There, cells detect 
the morphogen by receptors located on the cell surface and respond with patterns of gene 
expression that depend on the morphogen concentration. In this way, cell fate depends on 
the distance from the morphogen source and the interplay of several morphogen gradients 
originating from spatially distinct sources leads to the complex patterning of developing 
organisms. For this to be possible, the gene expression patterns must be formed with high 
precision and also need to be robust, i.e. they should not depend sensitively on parame- 
ters that are likely to fluctuate. Robust gene expression can result from the robustness of 
the morphogen concentration gradient p, Bl or can be achieved by other mechanisms j^. 
Interestingly, different species use very similar morphogens and a given species can use the 
same morphogen at different stages during development. An example are morphogens of 
the TGF-/3 super-family which exist in a wide range of organisms ranging from the fruit fly 
Drosophila to humans. 

The formation of morphogen gradients is still poorly understood. In particular, it is 
not clear, how morphogens are transported in the tissue adjacent to the source. Several 
possibilities have been proposed. Secreted morphogen molecules might simply diffuse in the 
extracellular space between cells and for some morphogens this indeed seems to be the 
dominant form of transport Q|. However, there is growing experimental evidence against 
simple diffusion for several morphogens jlll . 11^ . Other transport mechanisms, that have 
been proposed, include the bucket brigade, where molecules are "handed over" from one cell 
to adjacent cells 

Another possibility is suggested by recent experiments on the morphogen Dpp Dpp 
belongs to the TGF-/3 super family and patterns the developing wing discs of Drosophila. 
The wing disc is a sheet formed from one layer of cells that is the precursor of a wing. Fluo- 



2 



rescently labeled Dpp was tracked in wing discs where it is secreted by a narrow strip of cells 
from where it spreads into the adjacent tissue. If all cells adjacent to the secreting cells are 
defective in the internalization of the morphogen into the cell, Dpp gradients extend only 
about two cells from the source as compared to 30 cells in a non-mutant fly. If only cells in 
a small patch close to the source show this defect, a pronounced transient depletion of the 
morphogen is observed behind the mutant patch during gradient formation. This suggests 
that Dpp is not simply diffusing through the extracellular space. It has been suggested that 
these observations can be described by a scenario where Dpp diffuses in the extracellular 
space, binds to and is released from receptors [l^. However, an analysis which takes into 
account more recent exp erimental results indicates that free diffusion alone cannot account 
for the observations [l6[. An alternative suggestion is that long range transport of Dpp is 
generated by transcytosis, i.e. by repeated rounds of morphogen binding to cell surface re- 
ceptors, internalization into the cell and subsequent externalization and release of the ligand 
from the receptor at a different point on the cell surface A simple model to describe 
transcytosis is illustrated in Fig. ^ Ligands spread from a stripe of secreting cells into a 
two-dimensional tissue. However, in many situations a one-dimensional description of ligand 
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Transcytosis, can 



proflles as a function of the distance x to the source is appropriate 
be characterized by rates fcon and fcofr for receptor binding and unbinding of free ligands, as 
well as 6int and 6cxt for internalization and externalization of ligand-bound receptors. Inter- 
nalized morphogens are degraded with rate 6deg • The free diffusion of ligands in the narrow 
space between cells is characterized by the diffusion coefficient Dq. Diffusion alone does not 
lead to efficient long-range transport if the diffusion length = (-Do/cdeg)^^^ is smaller than 
one or several cell diameters. Here, Cdeg is the extracellular degradation rate. 

On scales large compared to the cell diameter a, this model together with the assumption 
that receptors are also synthesized and degraded in the cells, leads to effective transport 
equations for the number densities X{t,x) of morphogens and p{t,x) of receptors given by 

dtX = d,{D{X,p)d,X-Dp{X,p)d,p)-k{X,p)X (1) 

dtp = Z^syn(A, p) - Z/deg( A, p)p . (2) 

Here, D{X,p) is an effective diffusion coefficient and k{X, p) an effective degradation rate 
that depend on both the ligand and receptor concentrations. Furthermore, a term exists 
which describes ligand currents induced by gradients of the receptor concentration. We 



assume for simplicity that cells are non-polar such that externalization occurs with equal 
probability at any place on the cell surface. As a consequence, the resulting transport is 
non-directional. The kinetics of the receptor density is characterized by effective source and 
sink terms z^syn('^;P) and i^degi^^p)- The dependences of the effective diffusion coefficient 
D, the degradation rate k, the coefficients Dp, Usyn and z/deg on the ligand and receptor 
concentrations can be calculated explicitly In order to keep our discussion simple, we 
focus here on the case where the surface receptor number per cell R is maintained constant. 
Then, ligand transport is described by a single nonlinear diffusion equation 

dtX = d,{D{X)d,X)-k{X)X . (3) 



Explicit expressions of D{X) and A;(A) are given in [22|. The coefficients are displayed as 
functions of the ligand concentration A in Fig. |21 For small A, D{X) and k{X) approach 
finite values, and the transport equation becomes linear. For large A, -D(A) ^ Do + Ci/A^ 
with ci = abi^tkosr / {4:kon) , where r = R/a. The effective diffusion coefficient -D(A) decreases 
for large ligand concentration because receptors are increasingly occupied and not available 
for transport. A maximum of D occurs for intermediate values of A, if the ligand binding 
rate kon exceeds a critical value. This implies that for small A, D increases for increasing 
A. Similarly, the effective degradation rate exhibits for large A the asymptotic behavior 

k ~ Cdeg + C2/A where C2 = ^deg^int^^/^cxt- 

We solve Eq. Q for x > 0, representing the space adjacent to the ligand source, and with 
boundary conditions j{x = 0) = jo and j = for large x. Here, jo is proportional to the 
rate of ligand secretion in the source. In this situation, a graded ligand profile is generated 
and maintained in the steady state, see Fig. IHl 

The effect of the source is captured by the boundary condition on the current j = jo at 
X = 0, where j = —D{X)dxX, see Eq. (jH)). If initially at time t = ligand is absent A = 0, 
the ligand spreads for t > into the region x > 0, and builds up a gradient. In the steady 
state, the total ligand concentration decreases monotonically with increasing distance from 
the source. The steady state satisfies d^j + k{X)X = and is reached on a time scale of the 
order of 1/A;(A = 0). The corresponding ligand distribution X{x) obeys 

x=- f^'^ dX' D{X')/j{X') , (4) 
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where the hgand current j in the steady state is given by 

1/2 

(5) 

For small A, the ligand profile decays as A oc exp(— x/,^) with ^ = {D{0)/k{0)y^'^. In the 
case Do = 0, the steady state profile has a singularity at x = x*: 

A~ (x-x*)-i(-ln(x-x*))-i/2 (6) 

Choosing Jq > ensures x* < such that x* is not physically accessible. This singularity 
results from the asymptotic behavior of the current j for large A ^ A^: j(A)^ ~ H^tY + 
2edegCi ln(A/AT) + 2ciC2(l/Ar — 1/A). Here, At denotes a characteristic ligand concentration 
beyond which the limit of large ligand concentration holds. Therefore the current diverges 
as ~ 2ciedeg In A. Note, that for Cdcg = the current reaches for large A a finite maximal 
value jinax and the corresponding steady state ligand profile diverges as A ~ Ci/[jmax(a^ — a;*)]. 

The singular behavior of the steady state profile near x = x* has remarkable consequences 
for the robustness of gradient formation. Such a robustness can be quantified by an appro- 
priate response function of the system Qj. We define the robustness 7^(jo) := o-UodjoXi^X))^^ , 
where the function x{X) is given by Eq. (^. It becomes large if changes of jo have little 
infiuence on the ligand profile in the steady state. Here, a robustness of 7^ = 1 implies that 
under a 100% increase of jo the position at which the ligand profile attains any fixed value 
is displaced by about one cell diameter a. Thus for TZ>1, cells in the target tissue cannot 
detect variations of the ligand concentration due to significant changes of jo, while for TZ <1 
the steady state profile is strongly affected, see Fig. IHl Note, that TZ is independent of A. 
Indeed, TZ can be rewritten as 

7^ = a{jQdj^x)~^ = adxJo/D{Xo) = ak{Xo)\o/jo (7) 

where Ao = A(x = 0) and Eqs. (jlj) and © have been used. The robustness is thus completely 
determined by the ratio of the effective degradation rate and the ligand current at x = 0. 
High degradation rates and small currents lead to a robust gradient. Using the asymptotic 
behavior of the steady state profile for Z^o = 0, we find that the robustness increases rapidly 
for large currents jo as 7^ ~ jo^^e-'^/-'^ with j^ = l/2ciedog- For small jo, TZ — a/^ becomes 
constant. 
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The situation is different if Dq is finite. In this case, both k{X) and D{X) are nonvanishing 
and constant for large A and consequently Eq. Q becomes linear in this limit. The steady 
state profile for A ^ Ci/Dq behaves as A ~ exp(— x/^^). In the interval Ay < A < Ci/Dq the 
ligand profile is well described by Eq. (jHI), while for small A <^ At, A ~ exp(— x/^). As a 
consequence of finite Dq, the robustness approaches a finite value T^max = for large jo- 
For < a, this implies 7^ > 1 for large jo and cells in the target tissue are insensitive to 
variations of jo by a factor of two. The robustness is displayed in Fig. [Has a function of jo 
for different values of C,d/ci- 

Finally, we briefly discuss how the transport Eq. (jH)) is derived from a microscopic model. 
We denote the number of free ligands in the space between cell n and ri + 1 by L„ and the 
numbers of receptor bound ligands inside cell n by S*^*-* . Similarly, we introduce the numbers 
S^^ and 5*^'^ of receptor-bound ligand on the right and left surface of cell n, respectively, 
see Fig. [TJ The dynamics of these quantities is given by 

dtSl:^ = -KfiSl:^ + A;o„(| - - k,,tS^:^ + ^K^tS^^ (9) 

dtS^^ = -kosS^^ + Kn{^ - S'«)l„_i - 6i„t5i'^ + ^feextS-W (10) 

= -feext^i') + &i„t(5^) + Si[^) - 6dcg5«. (11) 

On large scales, we describe the ligand profiles by the densities l{x) = Ln/a, s{x) = S^''/a 
and s±{x) = [S'^'^±5'^''^]/a, where x = na. Because transport on large scales L is governed by 
time scales tl ~ L'^/D which are long compared to the relaxation time on the cellular scale 
Ta — a'^/D, three of the dynamic equations can be adiabatically eliminated. The effective 
Eq. (j2)) for the total ligand density A = / + s + S-|. follows by taking a continuum limit . A 
comparison of the full solutions of Eqs. (P)- (fTT|) to Eq. is shown in Fig. 13 demonstrating 
that the effective transport Eq. © captures the behavior of the microscopic model. 

In conclusion, we have developed a general theoretical framework for ligand transport 
and the generation of characteristic ligand profiles via transcytosis. Cell differentiation is 
triggered at threshold levels of morphogen concentration. The robustness of patterning has 
been demonstrated experimentally under conditions where Dpp was overexpressed |llL Il6, 
Q]. Reliable patterning can be achieved by robustly generating morphogen profiles. We 
have defined the robustness of steady state ligand profiles with respect to the rate of ligand 
secretion jo. We find that transcytosis naturally leads to a large robustness of the profile. 



Free extracellular diffusion reduces robustness significantly if the diffusion length exceeds 
several cell sizes. The origin of robustness in transcytosis are nonlinearities in the ligand 
current. A similar mechanism based on nonlinear degradation alone permits the robust 
formation of gradients in a free diffusion model U]. Robustness can also be achieved by 
other nonlinear effects ^J. 

In the general case where both ligand and receptor densities are taken into account, new 
phenomena appear Q]. In particular ligand currents can be driven by receptor gradients. 
Furthermore, ligand gradients are accompanied by a gradient in the receptor density as is 
indeed observed experimentally 11,|2^. In addition, if cells posess a polarity and define a 



direction in the tissue, transcytosis can lead to directed transport. Finally, it is straightfor- 
ward to generalize the concepts presented here to higher dimensions. 

The effects discussed here become important if the diffusion length C,d is smaller than the 
range ^ ~ 30a of the gradient. From experiments where endocytosis has been blocked in 
the target tissue, we estimate < 2a [11]. The rates of transcytosis have to be sufficiently 
large to guarantee gradient formation within the experimentally observed times of about 
A;(0)^^ ~8 hours. Assuming that the rates fcon-R, ^ofr, hnt, and 6ext are of the order of 1/min, 
we estimate D ~ a^/min. This diffusion coefficient is sufficient to generate a gradient over 
a length ^ ~ 30a during 8 hours. Numerical solutions to the transcytosis model (see Fig. 3) 
confirm this estimate. A rate of one per minute is larger than values measured in different 
but related systems j21] but is plausible for intracellular trafficking. 

The values of the rate constants are difficult to measure while the effective diffusion coef- 
ficient and the effective degradation rates are more easily accessible in fluorescence recovery 
experiments. Our work suggests that if transcytosis dominates transport, morphogen gra- 
dients are robust with respect to morphogen overexpression. This can be directly tested in 
future experiments. 
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FIG. 1: Schematic representation of transport by transcytosis in a chain of cells of diameter a 
indexed by n. The rates of ligand-receptor binding and unbinding, internalization and external- 
ization of ligand-receptor pairs are denoted kon, koS, hnt and 6ext- Degradation of ligand occurs 
inside the cells with rate 6deg and in the extracellular space with rate Sdeg- 
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FIG. 2: Effective diffusion coefficient D{X) as a function of ligand concentration A for Do/a^bdcg = 
10 (solid line) and Dq = (dashed line). Inset: effective degradation rate k{X) as a function of A for 
Sdeg = (dashed line) and edeg/&deg = 5 (solid line). Parameters are: 6int/^deg = ^ext/^deg = 3x 10^, 

konR/bdeg = 1-1 X 10^, kos/bdeg = 7 X 10^. 
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FIG. 3: Ligand densities A(x) in the presence of a source at x = at different times bdcgt = 
0.18,0.54,0.9,1.26 during gradient formation and in tlie steady state. Lines indicate solutions to 
Eq. ((31), while symbols indicate solutions to Eqns. ((5))- (|ll() for comparison. The robustness of the 
steady state profile is TZ !^ 470. Parameters as in Fig. |2lwith Dq = 0, jo/bdcgR = 70 and j = at 
x/a = 50. Inset: steady state ligand profiles for same parameters but jo/bdegR = 7 where TZ « 0.1. 
The profile (solid line) is strongly affected by halving (dotted line) or doubling (dashed line) the 
ligand current of the reference state. 
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FIG. 4: Robustness TZ of steady state ligand profiles as a function of the ligand current jo from 
the source for different values of the ratio of the diffusion length and the cell size a. Parameters 
as in Fig. [2 
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